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ABSTRACT 



Magnetic fluctuations in a non-magnetized gaseous plasma is revisited and 
calculated without approximations, based on the fluctuation-dissipation theorem. 
It is argued that the present results are qualitative and quantitative different form 
previous one based on the same theorem. In particular, it is shown that it is not 
correct that the spectral intensity does not vary sensitively with kcut- Also the 
simultaneous dependence of this intensity on the plasma and on the coUisional 
frequencies are discussed. 

Subject headings: plasma physics; magnetic fluctuations 



1. Introduction 



Fluctuations of physical quantities near zero frequency have been investigated by several 



authors since the papers of Johnson (1928) and Nyquist (1928). A general theory on the 



fluctuation-dissipation theorem, which will be the starting point of this paper, was developed 



by Kubo, R ( 1957). To the best of our knowledge, a concrete expression for the low-frequency 



spectrum of fluctuations of magnetic fields in a thermal plasma was obtained for the first time 



by Tajima, T.; Cable, S; Shibata, K. & Kulrsrud, R.M. (1992). They found a peak around 
u = magnetic fluctuation which was interpreted as the evanescent energy component of 
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electromagnetic fluctuations "screened" in a plasma below the plasma frequency. The impact 



of such result into the cosmic microwave background was then investigated by Tajima, T.; 



Cable, S; Shibata, K. & Kulrsrud, R.M (1992). Although in these two references the authors 



claim that the fluctuations were rigorously computed, several approximations were indeed 
made and they were not able to get an unique formula covering the low- and high-frequency 
spectrum. The aim of this paper is to reevaluate the derivation of the spectrum of magnetic 
fluctuations avoiding approximations in the low-frequency region and also in the transition 
between the low- and the high-frequency spectrum. Several different behaviors between our 
results and the previous one, mainly in the low-frequency part of the spectrum, are found 
and discussed. 



2. The first predictions 



The fluctuation-dissipation theorem developed by Kubo, R. (1957) is able to deal with 
the thermal fluctuations inside a plasma in or near thermal equilibrium. The expression 
for the magnetic field fluctuation in an homogeneous isotropic non-magnetized equilibrium 
plasma was obtained by Tajima, T.; Cable, S; Shibata, K. & Kulrsrud, R.M, (1992) looking 
at waves in such a plasma. In an electron-positron plasma, for example, the magnetic 
fluctuations in wave number and frequency space is given by 



57r 



X 



(a;2 + r]^)k^c* + 2a;2(^2 _ ^2 _ ^2)^2^2 + [(^2 _ ^2)2 + ^2^2]^2 
where kg is the Boltzmann constant and Up and 77 are, respectively, the plasma and the 



collisional frequencies. Here uj^ = + +, and 



UJp 



'-fm 



being n the particle density inside the plasma, Z is the atomic number of these constituents, 
e and m are respectively the charge and the mass of the electron, and 7 is the Lorentz factor 
given by 7 = (1 -|- kT/mc^). This result can be integrated in d^fc = Ank'^dk to get (the 
Fourier transform) 
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The integral of eq. ([T]) over wave numbers shows a high wave number divergence. Ac- 



cording to Tajima, T.; Cable, S; Shibata, K. & Kulrsrud, R.M. (1992), this is expected 



since the derivation is based on classical fluid equations of motion and the constant collision 
frequency rj is independent of k. However, they prefer to carry on their analyzes in the 
simpler phenomenological approach. In order to overcome the large k dependence, they flrst 
take the limit — )■ and then they integrate over k to inflnity, which corresponds to the 
vanishing cross section of collisions as A; — )■ oo. This is a very delicate point and we will turn 
back to this point in Section [3j For both the high frequency and high wave number limits 
the authors emphasized that the expression of eq. ([T]) has a substantial value only where 
— c^A;^ — ~ 0. The combined high-frequency and high wave-number limits were get by 
letting — )■ 0. The expression for the low-frequency spectrum was obtained by breaking up 
the k integral into two intervals, by introducing a cutoff value kcnt, with Xcut = kcntc/oope- In 
the integration from to kcut, V was kept flnite while in the integral from fccut it was used the 
approximation r/ — )■ 0. The expressions obtained for the high and low parts of the spectrum 
was, respectively: 



= Ls(^) [ —^i— edk + —-^ _ 2)3/2 



and 



Stt e(^"^p-/'=s^V - 1 ' V c / J {00'"^ + r]''^)x^ + ■ 



(3) 



27re('''^p=/*''s^)'^ — 1 V c / V 
where 6 is the Heaviside step function, rj' = rj/cope, co' = u/upe, and Up = Up/upe. 
Finally the zero frequency limit of the magnetic fluctuations is give by 

uj^o Stt vr2(e^^''='^ /^s^ - 1) \ c J 7]' Jq 



At this point the frequency spectral intensity was plotted for a temperature T = 10^" K 
by requiring that the value of fccut (or Xcut) provide a smooth behavior at the joint between 
the low-frequency spectrum and the black-body spectrum. The choice was fccut ~ i^pe/c or 
(a^cut ~ !)• The result for other values of the temperature were presented in another paper 
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by Tajima, T.; Cable, S; Shibata, K. & Kulrsrud, R.M. (1992). The main claims by these 
authors was that the intensity of the spectrum does not vary sensitively with fc^ut and that 
near u = the spectrum goes like w"^. Let us now show our results. 



3. General Result 

We have integrated eq. ([T| over k analytically, without any approximation, by the partial 
fractions technique. The exact result for the indefinite integral over the wave number is: 
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with the following definitions: £^ = c/ujp is the Debye length, x = u/up, y = k/k^ = kc/uj, 
V = r/T/wp, a = l-{x'^+ri')~\ b = {t]' /x){x^+r]')-\ R±{y) = y^±2a'y+r', a' = ^{r + a)/2, 
b' = v/(r - a)/2, r = {a^ + b^/^ sin(3e/2) = {ba' + ab')r-^/\ cos(30/2) = {aa' - 66')^"^^', 
and 
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where Xp = k^T/hvp. 



Note that the general result shows only a linear divergence in k restricted to the first 
term of eq. (|5|. This term, however, cannot be simply discarded by doing the limit r] ^ 
before we integrate over k, as did by Tajima, T.; Cable, S; Shibata, K. Kulrsrud, R.M. 
[| fl992| , smce it plays a very important role when the strict limit u/up — >■ (which we prefer 
to indicate from now on by w <^ Up) is to be considered, even when k is large. Indeed, if 
we discard it for all large values of k it can be shown that the limit x <^ 1 of S{uj) will be 
negative. Therefore, our result for the definite integral can be put in the form 



S{u) = 5*1 + 5*2 + 5*3 



where 



= Uujb / dy; S2 



{a'b + ab'); S3 



{aa — bb') 



(7) 
(8) 
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The term 5*1 will be taken as 5*1 = u^bi/cnt with ?/cut as large as we want. This will render 
the confrontation with the previous result easier. We note that for k > 10^^ cm~^, 5*2 — t- 0. 
These three terms can be written as functions of the original parameters, considering that 



and 



aa 



bb' 



iJ^ + if 



(a;2 + ^2) [(^2 _ ^2)2 + ^2^2] + ^(^2 + ^2 _ ^2) 



1/2 



+ ^2) [(^2 _ ^2)2 + ^2^2] _ ^(^2 + ^2 _ ^2) 



.1/2 



y2a;3/2(^2 + ^2)3/2 

Compared to eqs. ([2])-(|4]) it is immediately evident how our result is different from those of 
equations, showing a much more complicated dependence of the frequency spectrum on the 
variable a;, and on the parameters ujp and 77, which dependencies on plasma temperature are 
shown, respectively, in Figures [1} 




T(K} T(K) 



Fig. 1. — (a) Plasma frequency and (b) plasma collison frequency, both as a function of 
plasma temperature. 

Our full result S{uj) is shown in Fig. |2| where both the normalized (gray curve) and non 
normalized (green curve) spectral intensities are given, for kcut = 10^'' cm"^ e T = 10^° K. 

In this figure we have also plotted each one of the terms that contribute to S{u)) are 
shown in different colors. The deep we see in this figure near x = 1 tends to disappear as 
kcut goes to high values; as kcut decreases up to ~ 10^^ cm~\ the ordinate of the deep tends 
to zero. 

The detail of the non normalized spectral intensity near a; ~ is shown in Fig. [3} It 
shows naturally a smooth behavior between the low-frequency spectrum and the blackbody 
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10 '= 10" 10 '° lO' lO " lO ' 10 ' 10 ^ 10^ lO ' 10 = 10 ' 1 10 10= lO' 10^ X 



Fig. 2. — Non normalized Siio) and normalized (S'(a;)/S'o) spectral intensities showing the 
different components of the spectrum, S\, S2 and 6*3, as given by eq. ([s]). 



spectrum, which is constructed by hand in |Tajima, T.; Cable, S; Shibata, K. fc Kulrsrud, 



R.M, (1992). 



If we compare this graph to the correspondent one shown in Fig. 1(a) of the paper 



of Tajima, T.; Cable, S; Shibata, K. & Kulrsrud, R.M. (1992), we see that our result, for 
X > 1, is a classical blackbody radiation spectrum which goes to zero at x ~ 10^ while their 
blackbody spectrum has a much greater width (more than two orders of magnitude). 

Another important qualitative difference between the non approximate and the approx- 
imate results is that we found a very peculiar oscillations in S{u) for the x ^ 10~^ region of 
the spectrum, as can be seen from Figure |4j These oscillations occur in an x region where 
the classical blackbody spectrum still have a significant value; there is however, in this case, 
a strong interference in the total spectrum S{uj), eq. ([T]), due to a change of sign of the 
function 83(00). Such a kind of behavior was found just for cut-off values of the order of 
kcut = 10^° cm~^. For values of kcut greater than this one such fluctuations disappear. In 
any case, this feature confirm our statement that the result can vary sensitively with kcut, 



contrary to what was sustain by Tajima, T.; Cable, S; Shibata, K. & Kulrsrud, R.M (1992). 



Finally, we have studied the behavior of S{u,ri) by varying u e rj. The result is shown 



- 7- 




Fig. 4. — Detail of the spectral intensity fluctuations for kcut = 10 



iio ^^j^ i_ 
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in Figure |5} Notice that just the peak of the zero-frequency plasma spectrum depends on rj 
(the blackbody part remains unchanged). Indeed, the spectral intensity varies two orders of 
magnitude by varying t] by two orders too, namely, it goes from S{uj, rj) ~ 10®, for r] ~ 10~®, 
to ~ 10'', for 7] ~ 10~^. Thus, our result indicates that, when one goes backwards in time, 
temperature grows, dynamo action is enhanced (since r] goes down), and the resonance peak 
of the zero-frequency plasma peak goes down. 




Fig. 5. — Dependence of the spectral intensity on rj and oj. 



Discussions 



In this paper we have computed the spectrum of magnetic fluctuations of an homoge- 
neous cosmic plasma avoiding any approximations. Several different behaviors between our 
results and the previous one obtained by |Tajima, T.; Cable, S; Shibata, K. &: Kulrsrud, 



R.M. (1992), mainly in the low-frequency part of the spectrum, are found and discussed. 



It is important to stress that the exact result indicates that the peak of the zero-frequency 
spectrum can indeed vary sensitively with the cut-off value (kcut)- 



In the light of this new result, and following the papers of Tajima, T.; Cable, S; Shibata, 



K. & Kulrsrud, R.M. (1992) and Caruso, F. & Oguri, V. (2005), the problem of establishing 



an upper limit for fractal space dimensionality from COBE data can be revisited. 

Our results can still be improved towards cosmological applications by computing the 
Fourier transformed volume element d^k in terms of curved Riemannian space embedded 
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into general relativistic spacetime. This will allow us to address the problem of dynamo 
action in Einstein cosmology. The dynamo effect in plasmas is a competitive effect between 
convection of the cosmological fluid and the plasma resistivity. This is the reason why is so 
interesting to consider the relation between the dynamo action 7 and the plasma resistivity 
and its frequency. Recently, some new investigations on this directions were addressed by 



Rubashnyi, A.S. & Sokoloff, D.D. (2010), by using the magnetic field correlation tensor in 



space of negative curvature, and by de Souza, R. & Opher, R. (2010) in the case of positive 
curvature. Also recently, Garcia de Andrade, L.C (2010) has obtained a constraint on 



dynamo action from COBE data, using two-dimensional spatial sections of negative curvature 
of Friedmann universe, based on the general relativistic magnetohydrodynamic equation 



derived by Marklund, M. fc Clarkson, C.A. (2005) 



(1996) 



In near future, by using the formula of |Kulsrud, R.M; Cen, R.; Ostriker, J. P. fc Ryu, D 

deM 



Cm 



7 



where 7 is the magnetic field growth rate (dynamo action), we expect to compute 7 in terms 



of the magnetic plasma resistivity (77) and plasma frequency 
estimate the variation of dynamo action in terms of r] and Up. 



In this way we shall 
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